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ABSTRACT 

Closed form and finite element models are developed 
to predict the device response of aluminum electrostatic 
torsion micro-mirrors fabricated on silicon substrate for 
space applications at operating temperatures of 30K. 
Initially, closed form expressions for electrostatic pres- 
sure and mechanical restoring torque are used to pre- 
dict the pull-in and release voltages at room temper- 
ature. Subsequently, a detailed mechanical finite ele- 
ment model is developed to predict stresses and ver- 
tical beam deflection induced by the electrostatic and 
thermal loads. An incremental and iterative solution 
method is used in conjunction with the non-linear finite 
element model and closed form electrostatic equations 
to solve the coupled electro-thermo- mechanical problem. 
The simulation results are compared with experimental 
measurements at room temperature of fabricated micro- 
mirror devices. 

Keywords : micro- mirrors, simulation, cryogenic, elec- 
trostatic, thermal. 

1 INTRODUCTION 

NASA Goddard Space Flight Center (GSFC) is de- 
veloping two-dimensional tilt-able aluminum Micro- Mir- 
ror Arrays (MMA) for Earth and space science applica- 
tions, such as actively controlled masks for spectrom- 
eters. Due to science requirements, each micro-mirror 
element must be a factor of six larger than the Dig- 
ital Micro-mirror Devices (DMD) developed by Texas 
Instruments, while maintaining similar rotation angles, 
resulting in much larger gap distances for electrostatic 
actuation. In addition, these devices must operate at 
30K, resulting in potentially severe thermal mismatch 
stresses. 

Our objective is to enable optimal design of the GSFC 
MMA through detailed modeling and simulation of a 
single micro-mirror element shown in Fig. 1. In this 
work we present improved theoretical and finite element 
models to simulate the device response at ambient and 
cryogenic temperatures. In particular, a theoretical mod- 
el similar to that presented in Ref. [1] is developed, but 
extended to account for large angles of rotation. Also, 
a finite element model is developed using shells instead 



Figure 1: Single micro- mirror element. 

of bricks for high aspect ratio members to more accu- 
rately account for the vertical deflection and stiffness of 
the torsion bars. The analytical and finite element re- 
sults are compared with results from fabricated micro- 
mirrors. 

2 ELECTRO-MECHANICAL 
SIMULATION 

In this section we present approximate equations gov- 
erning the electro- mechanical response of a micro- mirror 
element. These equations are developed neglecting the 
vertical deflection of the torsion bars. As a result, we 
also present a finite element model, which can be used 
to more accurately estimate local stresses. 

The models are developed for the set of parameters 
shown in Fig. 2. In this figure the coordinate system is 
aligned with the torsion bars as shown. The torsion bar 
dimensions are length Lb ar , width Wbar, and thickness 
tbar . Three key points are labeled as a, 6, and c. Their 
respective x locations are x a , ^6, and x c . The only y 
location needed in the simulation is Under an applied 
voltage the yoke rotates through angle 0 about the y 
axis. 



Figure 2: Post, torsion bar, and yoke geometry and pa- 
rameters. 

2.1 Analytical Model 



The micro-mirror motion is driven by applying an 
electrostatic field between two plates offset from the tor- 
sion bar central axis. The two plates are non-parallel 
in general, and the electrostatic pressure distribution is 
non-uniform. The resulting pressure distribution, de- 
rived using the Laplace (Ref. [2] pg. 34) and Maxwell 
Stress Tensor (Ref. [2] pg. 261) equations, is 
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Figure 3: Electrostatic and mechanical torque versus 
yoke inclination angle. 

where 0 is the twisting angle and corresponds with the 
relative inclination angle 0 between the yoke and elec- 
trode pad. The term G is the shear modulus, and Lbar 
is the length of one of the torsion beams. The factor 
of 2 accounts for both torsion beams. Finally, K is the 
torsional stiffness, which is [3] 


where £o is the permittivity of air, 6 is the inclination 
angle, and p(x,0), for the coordinate system shown in 
Fig. 2 is 


In the above equations, V is the applied voltage, do is 
the initial distance between the two electrostatic plates, 
and x is the coordinate relative to the system shown 
in Fig. 2. These equations neglec t electrostatic edge 
effects, but are not limited to small inclination angles. 
The total electrostatic torque is 
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The above equation is for a beam with solid rectangu- 
lar section, of width 2a and thickness 2b. In this case, 
a = iVbar/2, and b = t^ aT / 2. Equations (3) and (5) are 
plotted in Fig. 3, for L^ ar = 47.0 /zm, wtar — 6.0 pm, 
tbar = 0.3 pm, x a = 12.7 pm, x.b = 29.5 pm, x c = 
35.5 pm, yi } = 34.5 pm, and do = 5.8 pm. These pa- 
rameters yield a theoretical contact angle of 10°. 

At static equilibrium, the electrostatic and mechan- 
ical torques are coupled through the equation 
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where the edges of the yoke between points a and b are 
approximated by the line y(z) = Xb + yb — x. Equa- 
tion (3) can be be expanded to 
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The reacting mechanical torque can be computed us- 
ing beam equations. In this case the mechanical torque 


is 


T m (6) = 2 
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Tm(0) = T e (V,d). (7) 

Points where the curves intersect represent solutions to 
equation (7). As shown in Fig. 3, there are three types of 
intersecting points, labeled A, B, and C . For a low equi- 
librium voltage V e , the curves will intersect in the region 
near A. As the voltage increases the intersection moves 
up the dashed curve. When the voltage exceeds the crit- 
ical snap-on voltage V s at point B y the yoke unstably 
rotates until it contacts the substrate. The yoke is held 
in the contact position by a voltage much greater than 
needed for static equilibrium in the region near point 
C . As a result, the yoke will remain in contact until the 
voltage falls below the equilibrium release voltage V r for 
a given contact angle in the region near C. At this point 
the voltage will be too low to maintain equilibrium in 
region C and the yoke will snap back to region A. 






Figure 4: Solutions to equation (8) for various torsion 
bar widths. 

Equation (7) can be used to determine V in terms of 
6 as follows 
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The above equation is plotted in Fig. 4 for various 
bar widths uw, and the points A, B , and C are labeled 
consistent with Fig. 3. Equation (8) can be used to 
determine V e and V r - At the snap-on voltage, both V s 
and 6 S are unknown. However, these values are located 
at the local maximum. For the dimensions used in Fig. 3 
the resulting snap-on voltage and angle of inclination are 
V s = 24 V and 6 S — 4.8°, respectively. For a contact 
inclination angle of 6 r = 10°, the corresponding release 
voltage is V T = 12 V . 

Maximized solutions to equation (8) can also be used 
to study the sensitivity of the snap-on voltage V s to 
various geometric parameters. Fig. 5 show's the variation 
in V s for four key geometric parameters. This figure 
shows that the initial plate separation r/o has the most 
impact on V s . Increasing the bar width w^ar has nearly 
the opposite effect as increasing Xf,. The parameter x a 
has much less effect than £5, which matches intuition. 


2.2 Finite Element Model 

In this section, we present results from a detailed me- 
chanical finite element model of the key structural com- 
ponents in conjunction with equation (1). This model 
more accurately represents the system stiffness, and ac- 
counts for the vertical deflection neglected by the ana- 
lytical model. 

The mechanical finite element model consists of the 
posts, torsion beams, and yoke as shown in Fig. 6. Brick 
elements are known to significantly over estimate the 
stiffness of high aspect ratio structures. As a result, the 


Figure 5: Sensitivity of the snap-on voltage (K s ) to var- 
ious geometric parameters. 

torsion beams and yoke are modeled with shell elements. 
The mechanical properties of bulk aluminum are used 
for the post, beam, and yoke. The base of the posts are 
fixed, and equation (1) is used to compute the pressure 
distribution for a given voltage and angle of inclination. 

The pressure distribution is a function of the orien- 
tation of the yoke. As a result there is a highly non- 
linear and coupled interaction between the electrostatic 
and mechanical behavior, which we simulate with a nu- 
merical procedure. In this procedure, the voltage is ap- 
plied in small increments, and the solution iterates be- 
tween the electrostatic and mechanical solutions until 
achieving convergence. During the voltage increment 
that spans the snap-on voltage V s , the solution for 6 in- 
creases unstably until the yoke contacts the substrate. 
Consequently, the mechanical finite element model must 
account for contact between the edge of the yoke and the 
substrate. The above procedure was implemented with 
a program written in the PERL computer language [4], 
which calls UAI/NASTRAN [5] as needed to solve the 
finite element equations. 

The finite element (FE) and analytical models were 
used to plot 6 versus input voltage in Fig. 7 for Wbar = 4 
and 6 fim. Similar micro-mirrors were fabricated and 
tested [6], and the measured snap-011 voltage V s of ap- 
proximately 14V for iVbar “ 6 l irn is a ^ so included. The 
FE results show that the contact angle is less than 10° 
due to vertical deflection of the torsion beams. The FE 
results over-predict the measured snap-011 voltage, and 
the analytical results over-predict the FE results. The 
first discrepancy may be due to a mismatch between the 
model and experimental initial offset do, wdiich is diffi- 
cult to measure experimentally, and the finite element 
model may be too coarse. The second discrepancy is 
due to the vertical deflection, which is not included in 
the analytical model. Despite these discrepancies, the 






Figure 6: Finite element mesh of the yoke, torsion bar, 
and post assembly. 

analytical and finite element models predict the same 
trends, which validates the slopes presented in Fig. 5. 

3 THERMAL SIMULATION 

While the fabrication temperatures of the micro-mir- 
rors are 20 °C or above, the operating temperature is 
30K. The resulting change in temperature of -2G3K, cou- 
pled with the Coefficient of Thermal Expansion (CTE) 
mismatch between the aluminum structure and silicon 
substrate results in significant internal stresses. 

The average Young’s moduli of aluminum and silicon 
from 293 K to 30# are approximately 71.7 x 10 9 GPa 
and 175 x 10 9 GPa , respectively. While the thickness of 
the torsion bars are approximately 0.3 fim, the substrate 
thickness is on the order of 600 fim. Hence, the substrate 
will contract independent of the supported structure. 
The resulting net thermally induced relative displace- 
ment between the two posts is 

& net — (^aZ (^) 

where a is the CTE averaged over the specified temper- 
ature range, and Lt oia i is the total structure length. In 
this case n n i = 19.3e— 6 ppm/ A", a sl . — 0.82c— 6 ppm/ K, 
and L total = 135 jim. The resulting relative displace- 
ment is S ne t — 0.656 fim. 

The above lateral displacement results in an out-of- 
plane bending moment on the torsion bars and non- 
linear geometric behavior due to the offset between the 
yoke and torsion beam mid- planes. As a result, 
was applied to the base of one post in the finite ele- 
ment model while the other post remained fixed, and 
a non-linear geometry finite element analysis was com- 
pleted using UAI/NASTRAN. The solution predicts a 
peak Von Mises stress near the ends of the torsion bars 
of approximately 445 MPa due to the combined axial 
tension and bending. This value exceeds the room tem- 
perature yield strength of aluminum by 35%. However, 
fabricated micro-mirrors have survived cryogenic test- 
ing, indicating that local yielding in the torsion bars is 
likely. 



Figure 7: Voltage versus yoke inclination for two differ- 
ent torsion bar widths. 

4 CONCLUSION 

Analytical and numerical models of aluminum micro- 
mirrors subjected to cryogenic temperatures were pre- 
sented. While the model predictions are in coarse agree- 
ment with measured values, the analytical and numeri- 
cal models predict similar trends due to geometric vari- 
ations. The differences in model and measured data fall 
within expected shortcomings due to modeling assump- 
tions, and can be corrected with further model develop- 
ment and correlation. In particular the analytical model 
can be extended to account for the vertical deflection. 
The finite element mesh can be refined and the initial 
plate separation do can be adjusted to correlate with 
experimental results. In future work, the thermal and 
electro- mechanical models can be merged to predict the 
response at cryogenic temperatures. 
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